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Predictions for the difference of proton-antiproton and proton-proton 
cross sections in the heavy quark production within LO DGLAP analysis 
together with the ln 2 x terms resummation are presented. An important 
role of the double logarithmic ln 2 x corrections in a case of the large rapidity 
gap between a quark and an antiquark is discussed. 
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1. Introduction 

Tevatron, LHC and SSC experiments enable to investigate a new kine- 
matical region of hadron collisions, where the hadron center-of-mass energy 
y/s is much larger than the parton momentum transfer Q. Large enough 
scale Q 2 justifies the use of perturbative QCD (PQCD) and therefore hadron 
scattering data could be a test not only of the Regge theory (soft physics) 
but of different PQCD approaches e.g. DGLAP or BFKL too. The high 
energy hadronic processes depending upon the values of the momentum 
transfer squared Q 2 can be divided on the soft (with small Q 2 and small x 
involved) and the hard ones (at large Q 2 > lGeV 2 ). The high energy be- 
haviour of the total hadronic cross sections in the Regge theory is described 
by two contributions: one is close s ' 1 , the other is close s -0,5 . The first term 
~ s results from soft pomeron exchange while the second one ~ s -0,5 is 
identified as corresponding to p, u, f'2 or ai exchange. Thus the Regge poles 
theory predicts a decrease of a tot (s) as s — > 00. Experimental data show 
however that at s ~ lOOGeV 2 the total hadronic cross sections have rather 
a weak energy dependence i.e. they slowly (logarithmically ~ ln 2 s) increase 
with energy at larger energies. On the other side PQCD predicts very rapid 
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increase of parton distribution functions (PDFs) at small x (i.e. large s) 
region. In the leading ln(l/x) BFKL approach one obtains bare Pomeron, 
which causes that cross sections increase as a power of energy, what is the 
evident unitarity violation. Very interesting both experimentally and the- 
oretically are so called semi-hard processes (SHP), which involve relatively 
large transverse momenta Q 2 (Q 2 » A 2 ~ 0.04GeV 2 ) and soft parton colli- 
sions with low momentum fraction x (x < 10~ 2 ) as well. In this way SHP 
open the description area in which the Regge theory and PQCD adjoin. Un- 
fortunately, in this semihard regime there is a uncertainty in determination 
of the small x dependence of parton distributions. In order to avoid this 
uncertainty in predictions of jet production cross sections it is necessary to 
measure the two-jet inclusive cross section in hadron collisions with a rapid- 
ity gap between these two jets. In such processes the dominant part is an 
exchange of BFKL Pomeron between partons. Events with a rapidity gap 
i.e. a region of rapidity in which no particles are found occur if at the parton 
level the colour singlet (qq or gg) is exchanged in the t-channel. Though the 
Pomeron (gg) exchange mechanism is the leading one, the second contri- 
bution (qq) coming from the mesonic reggeons is interesting too. Analysis 
of the mesonic reggeon exchange enable to test double logarithmic effects 
in PQCD. Elimination of the background origing from Pomeron, reggeized 
gluon and odderon in QCD analysis of hadron collisions gives possibility to 
determine the appropriate cross sections via valence quark- antiquark elas- 
tic scattering terms. The difference of the cross sections of PP and PP 
collisions incorporating ln 2 x approximation in valence quark distribution 
functions is not too small (> 100 pb) and can be measured e.g. in Tevatron 
and LHC. 

In our paper we show predictions for the heavy quark production cross 
section in hadron collisions. Taking into account the difference a(PP — ► 
hh) — a(PP — > hh) we are able to investigate the pure mesonic exchange 
picture modified by PQCD ln 2 x effects. In the next section we shall briefly 
recall the cross section formula for the production of heavy quarks in hadron- 
hadron collisions. We shall show how to eliminate all contributions except 
one origing from the mesonic reggeon exchange. Then in point 3 we shall 
present PQCD modification of the Regge model approach in which the ln 2 x 
effects in valence quark distributions are taken into account. We calculate 
cross sections a(PP — > hh) — a(PP — ► hh) and da(PP — > hh)/dAy — 
da(PP — > hh)/dAy for heavy quarks charm and bottom production using 
GRV input parametrizations of valence quarks and including the ln 2 x effects 
into the DGLAP evolution of these distribution functions. All our numerical 
results we shall present in section 4. We shall emphasize the role of the 
rapidity gap in our approach. Finally in conclusions we shall roughly discuss 
future experimental hopes for observation of ln 2 x effects in hadron-hadron 
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collisions. 



2. Production of heavy quarks in hadron-hadron collisions. 
Valence quark- ant iquark elastic scattering as the only 
contribution to a(PP — > hh) — a(PP — > hh) cross section. 

The total cross section for heavy quark production in hadron-hadron 
collisions via factorization theorem can be written in the form 1 : 

a(AB -> hh) = J dxiPi/A(xi,Q 2 ) J dx 2 Pj/ B (x2,Q 2 ) J ~^di 

4m? 4m? ^min 
. h h 



where A,B denote colliding hadrons, h(h) is a heavy quark (antiquark) 
with mass nih, i,j run over all partons i.e. quarks and gluons: (i,j) = 
(9,9), (q a ,q a ), (q a ,q a ),a = u,d,s. Pi/ A (x,Q 2 ) is the parton of i sort (quark, 
antiquark or gluon) distribution in hadron A, are the longitudinal mo- 
menta fraction of partons with respect to their parent hadrons. Q 2 is the 
parton momentum transfer squared (Q 2 ~ 4m? ). da/dt is the differential 
partonic cross section, which in a case of quark-antiquark subprocess in LO 
approximation has a form 2 : 

msi^m _ [(i _ mtf + ( , _ mlf + 2<i] (2 . 2) 

a s (Q 2 ) is the strong coupling constant; s,t,u are the well known internal 
Mandelstam kinematic variables and 

s + t + u = 2m 2 h (2.3) 

Hence (|2.2|) takes a form: 

da{qq^hh) _ ±ixa 2 s {Q 2 ) 4 „ 2 „ 



eft 9s 4 



(2m£ + s 2 + 2t J + 2st - 4m£f) (2.4) 



After integration with respect to t one can obtain the total partonic cross 
section, which in a case qq part has a form: 

J (if 

The lower and upper limits in the above integral are: 

t min = ^-(l + P) 2 (2.6) 
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t™„ = ^(i-/?r (2.7) 

with 

P 2 = 1 _ M (2 8) 

S 

and hence one can simply obtain: 

Q 2 ) = l! ^ 1 /3(3 - /3 2 ) (2.9) 

The total cross section for heavy quark production is of course dominated 
by the gluon-gluon gg — > hh subprocess and if we want to deal only with 
a part coming from the quark-antiquark scattering we should consider suit- 
able difference of cross sections. As it has been exactly described in [3], the 
simplest way to eliminate the strong background due to the exchange of 
Pomeron, reggeized gluon and oderon is to regard PP (proton-antiproton) 
and PP (proton-proton) collisions. All these mentioned above backgrounds 
cancel if one takes the difference of the cross sections a{P P) — a{P P) . Then 
the difference of the PP and PP cross sections for the heavy quark pro- 
duction is purely represented by valence quark-antiquark elastic scattering 
(with mesonic reggeon exchange in the Regge theory). One can write the 
differential cross section as [Sj: 

da(PP -> hh) da(PP -* hh) _ 1 ^ Q 2 )Aq (x Q 2 ) 

dx\dx2di dx\dx2di 2 ~ " P ^ P 2 ' 

+q f/P (x 1 ,Q 2 )Aq f/P (x 2 , Q 2 )\~^ 



with 



(2.10) 

Aq f/P (x, Q 2 ) = q f/ p{x, Q 2 ) - q f/P (x, Q 2 ) (2.11) 
Aq f/P (x, Q 2 ) = q f/P (x, Q 2 ) - q f/P (x, Q 2 ) (2.12) 

By summation over flavours / in parent hadrons (P or P) in (|2.1U|) one 

gets: 

da(PP -> hh) da(PP -> hh) _ 1 ^ ^ ^ 2 _ q2\^_ 

dx\dx2di dx\dx2di 2 ^ P 2 ' di 

(2.13) 

Because 

Qval/pix, Q 2 ) = q va i/ P (x, Q 2 ) (2.14) 
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one has: 



da(PP -> hh) da(PP -> hh) _ 1 2 2 



dx\dx 2 dt dx\dx 2 dt 2 



■^[u va i{xi,Q )u val (x 2 ,Q 



+d val (x 1 ,Q 2 )d val (x 2 ,Q 2 )}^ (2.15) 

at 



or after integration with respect to t: 

da(PP -» hh) da(PP^hh) 2-kol 2 s {Q 2 ) 



m - P 2 



dx\dx 2 dx\dx 2 27s 

x[u va i(xi,Q 2 )u va i(x 2 ,Q 2 ) + d va i(xi,Q 2 )d va i(x 2 ,Q 2 )} (2.16) 

where ft 2 is defined by ()2.8j) and 

Q 2 « 4m£ (2.17) 
Then the difference of the total cross sections is: 



[a{PP -> M) - a(PP -> fefe)](s) = } / dxi / dx 2 

4m? 4m? 
ft 

x K a ;(xi, Q 2 )u va i(x 2 , Q 2 ) + d va i(x 1 ,Q 2 )d val (x 2 , Q 2 )] 

(2.18) 

with internal kinematical variable s: 

s = sx\x 2 (2-19) 

In the formula Q2.18JI we have contributions from all possible rapidity gaps 
Ay between a quark and an antiquark: 

Ay = y(h) - y(h) « ln-^ (2.20) 

and hence 

< Ay < In-^r (2.21) 
~ 4ml 

In a case of Tevatron with y/s = 1.8TeV and for cc production with m c = 
1.43GeV theoretically 

< Ay < 13 (2.22) 
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However the interesting double logarithmic ln 2 x effects are important when 
the rapidity gap Ay is large: 



(Ay) z = In 2 tt ~ ln 2 x (2.23) 

Am 2 



To see this one should take into consideration the differential cross section: 

da(PP^hh) da{PP -> tih) _2^a 2 s {Q 2 ) f dx x . 2 

— p(o — p ) 



dAy dAy 27s J x\ 



4m? e As 



x [u„ a j (si , Q 2 )u m/ ( — * , Q 2 ) + d„ o/ (xi , Q 2 )d val ( — * , Q 2 )} 

s s 

(2.24) 

with (3 function of Ay: 

(3 2 = 1 - e"^ (2.25) 
From (|2.25|) one can simply notice that for a large Ay /3 2 is close 1.0: 

(3 2 = 1 as Ay — ► oo (2.26) 
and hence (|2.24|) for a large Ay can be expressed as: 

da(PP^hh) da(PP^hh) 4Tra 2 {Q 2 ) } dxi 



dAy dAy 27s J x\ 



4m? e A V 
h 



a 2 Ay 4 2 Ay 

x[tw(*l,Q 2 )«™*(— ,Qp) + d vol (x 1 ,Q 2 )d wd (—* iQ 2 )} (2-27) 

s s 

According to (|2.23() . with larger rapidity gap Ay the ln 2 x effects in the 
a(PP — > hh) — a{PP — * hh) would be better visible. These effects we 
introduce via suitable modification in DGLAP LO evolution equations for 
valence quarks u va i and d va [. The perturbative QCD approach with an 
appearance of the double logarithmic terms ln 2 x can modify the Regge 
model. We discuss this in the next section. 



3. Double logarithmic effects ln 2 x for flavour nonsiglet parton 

distribution functions. 

The high energy s limit corresponds by definition to the Regge limit. 
In the Regge pole model jl] an amplitude of high energy elastic scattering 
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T(s,t) ~ s a ^\ where a(t) is the Regge trajectory, which is equal to spin J 
of the corresponding particle at t = M 2 (M is the mass of the exchanged 
particle). Hence the total interaction cross section, which by the optical 
theorem is connected to ImT(s, 0) can be written as a sum of the Regge 
poles contributions: 

a tot (s) = Y,h(0)s ah(0) - 1 (3.1) 

k 

where «fc(0) are the intercepts of the Regge trajectories and bk denote the 
couplings. The high energy behaviour of the total hadronic cross section 
can be described by two parts: first contribution is so called soft pomeron 
with intercept ~ 1.08 and the second one are the leading meson Regge tra- 
jectories with intercept or(0) rj 0.5. These reggeons correspond to p,u,f 
or a2 mesons exchanges. The flavour nonsinglet part e.g. valence quark 
distribution functions or F^ 3 = i 7 ^ — F^ is governed at small x (i.e. at 
high energy s = Q 2 (l/x — 1)) by 02 reggeon: 

Regge: q val {x) ~ x ~ a ^ (3.2) 

with a a2 (0) ~ 0.5. This behaviour is stable against the leading order 
DGLAP QCD evolution 5]. In other words the PQCD Q 2 evolution be- 
haviour of q va i at high energy is screened by the leading Regge contribution 
()3.2|) . But there is a novel effect concerning the flavour nonsinglet functions 
which is an appearance of the double logarithmic terms ln 2 x [H], [7j. For 
not too small values of the QCD coupling a s this contribution is approxi- 
mately the same (or even greater) in comparison to the contribution of the 
02 Regge pole: 

ln 2 x : q va i(x) ~ x _aJ (3.3) 



and 
where 



q = 2^ s (3.4) 

* s = — (3.5) 
At a typical value of a s (Q 2 = 4m^) for heavy quark production: 

a s {Q 2 = 4ml) = °- 30 ( 3 - 6 ) 

a s (Q 2 = 4m 2 b ) = 0.21 (3.7) 
u) becomes very close to a a2 : 

uj(a s (Q 2 = 4m 2 c ) » 0.5 (3.8) 

uj(a s (Q 2 = 4m 2 b ) w 0.4 (3.9) 



It must be however emphasized that even in the case of the PQCD sin- 
gularity generated by the double logarithmic ln 2 x resummation, this ln 2 x 
behaviour of flavour nonsinglet parton distributions is usually hidden be- 
hind leading ci2 reggeon exchange contribution. Only in the case of large 
rapidity gap Ay, where the intercept of the Regge trajectory a a2 (0) be- 
comes close to or even negative, it is possible to separate from the Regge 
background the ln 2 x behaviour of the considered cross sections. In our 
analysis we calculate the difference of cross sections (|2.18p and Q2.24JI using 
u va i(%,Q 2 ) and d va i(x,Q 2 ) obtained from DGLAP LO evolution equations, 
incorporating the double logarithmic effects ln 2 x as well. This method 
combining DGLAP LO and ln 2 x approaches in the case of polarized flavour 
nonsinglet functions was presented in jHj, JO]- For unpolarized ones, 
which we now investigate, the suitable equations are the same, only the 
input parametrizations q va i(x,kQ) are different. In our approach we expect 
that the cross sections (|2.18|) and particularly this for the large rapidity gap 
(|2.24j) are governed by ln 2 x terms i.e. by the valence quark ladder, shown 
in Fig.l. The equation taking into account both DGLAP LO evolution and 
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Fig.l. A ladder diagram generating double logarithmic In (1/x) terms in the 
flavour nonsinglet part of heavy quark production cross section with large rapidity 
gap Ay. 
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ln 2 x effects for q va i{x,Q 2 ) functions has a form j^j: 



1 k 2 /z 



f(x,k 2 ) = /(°)(x,fc 2 ) + ^i J ^f(- z ,k> 2 ) 

X u2 

a s (k 2 ) k fdk /2 A } dz(z + z 2 )f(x/z,k' 2 )-2zf(x,k' 2 ) 



+ 



2n J k' 2 l 3Jz 1 - :. 

1 



t-- 



where 



h q + lln(l-x))f(x,k' 2 )] (3.10) 



f W( x ,k 2 ) = " s(fc2) [ 4 ] ^l±^M z l^ 2z ^M 



2tt 3 J z 1 — z 



x 



+ (i + !ln(l (3.11) 

The unintegrated distribution / in the equation (|3.1Ujl are related to the 
qval(x,Q 2 ) via 

Q 2 (V*-i) , v 
^(x,Q 2 ) = ^°J(x) + | (3.12) 



re 



where 



A;, , 



2 

«S0O = fc o) = / ^-/(x, A; 2 ) (3.13) 

o 

In above equations q va i(x,Q 2 ) denotes u va i(x,Q 2 ) and d va i(x,Q 2 ) distri- 
butions as well. In our calculations we use LO fitted GRV input 
parametrizations q va i(x, fog): 

A; 2 = lGeV 2 , A(n/ = 4) = 232MeV (3.14) 

G# : Uoaj(a;,jfeg = 1 GeV 2 ) = 2.872x~ a427 

x (l- 0.583x 0175 + 1.723x + 3.435x 3/2 )(l-x) 3 ' 486 (3.15) 

GRV : d val (x,k% = 1 GeV 2 ) = 0.448x" - 624 

x (l + 1.195x a529 + 6.164x + 2.726x 3/2 )(l-x) 4 - 215 (3.16) 

The numerical results for ()2.18j) and ()2.24|) cross sections we present in the 
next point. 



10 



4. Numerical predictions for cr(PP — > hh) — cr(PP — >■ hh)(s) and 
da{P dl7 hh) - da(P dl7 hh) incorporating DGLAP LO evolution 
and the ln 2 x effects in q va i functions. 

Our numerical results based on GRV (|3.14j) - (j3.1fi[) input parametriza- 
tions are presented in Figs. 2-7 and in Table I. In Fig.2 we plot xq va i input 
parametrizations for u and d valence quarks in proton at k 2 = lGeV 2 to- 
gether with their LO DGLAP + ln 2 x predictions at Q 2 = lOGeV 2 . Figs.3,4 
show the difference of total cross sections Q2.18JI which incorporates ln 2 x 
terms in q va i(x,Q 2 ) functions. Our calculations concern the charm quark 
(Fig. 3) and the bottom quark (Fig. 4) production. In Fig. 5 we compare LO 
DGLAP + ln 2 x prediction for the cross section (|2.24j) with pure LO DGLAP 
one at large Ay = 4.0. In Figs. 6, 7 we plot the difference of differential cross 
sections (j2~!?i|) with LO DGLAP + In 2 

x Qval( x i Q 2 ) f° r both charm and 
bottom production. The plots in Figs. 6, 7 are performed for different large 
rapidity gap Ay > 2.0. All plots in Figs. 3-7 are presented as a function of 
the center mass energy y/s. Finally in Table I we present the differential 
cross section 

da{PP -> hh) da{PP ->hh),_ ^a 2 s (Q 2 ) 

X\X 2 [ 



dx\dx 2 dx\dx2 27s 

x [u va i(x 1 ,Q 2 )u va i(x 2 , Q 2 ) + d va[ (xi, Q 2 )d va i(x 2 , Q 2 )] (4.1) 

for large rapidity gaps at different yfs and suitable x\ = x 2 '- 

2mh 4? (A <^ 

xi = x 2 = — -j^e 2 (4.2) 



From Figs. 3-4 one can read the shape of the cross sections difference 
a(PP — > hh) — a(PP — > hh)(s) as a function of y/s. An initial increase of 
its value with the increasing value of yfs results from the energy threshold 
for heavy quark production: 

s>s> Am 2 h ( 4 -3) 

where s is defined in (|2.19j) . When the energy becomes larger (y/s > 50GeV 
for charm and \/s > 150GeV for bottom production) because of a suppres- 
sion factor in the cross section formula (|2.18|) . a(PP) — a(PP) decreases 
with the further increasing of energy yfs. More interesting is analysis of the 

differential cross section ^^^1/^ ~ ^^dAy^ ■ When the rapidity gap Ay 
between a quark and an antiquark defined in (|2.20|) is large, according to 
()2.23j) high order corrections ~ (a s ln 2 x) n are important. It is shown in Fig. 5 
where we compare the LO DGLAP + ln 2 x predictions for the cross section 
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Fig. 2. Valence quark distribution functions xu va i and xd va i at input scale fc 2 , = 
lGcV 2 - dashed and at Q 2 = lOGcV 2 - solid. 
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Fig. 3. Predictions for a(PP -> hh) - a(PP -> hh)(s) incorporating LO DGLAP 
+ ln 2 x q va i{x,Q 2 ) for the charm quark production. 
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Fig. 4. Predictions for a(PP -> hh) - a{PP -> Wi)(s) incorporating LO DGLAP 
+ /n 2 a; q va i{x, Q 2 ) for the bottom quark production. 
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Fig. 5. Comparison of pure LO DGLAP predictions for dZ^[&(PP — * hh) — 
a(PP -> hh)] - dashed with LO DGLAP + ln 2 x ones - solid at Ay = 4.0 for 
the charm quark production. 
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X\ = x 2 


V(s)[GeV] 


Ay 


F [pb] 


0.54 


39 


4.0 


3776 


0.034 


630 


4.0 


3811 


0.091 


630 


6.0 


1090 


0.25 


630 


8.0 


236 


0.012 


1800 


4.0 


1628 


0.032 


1800 


6.0 


496 


0.087 


1800 


8.0 


143 


0.0015 


14000 


4.0 


277 


0.0041 


14000 


6.0 


90 


0.011 


14000 


8.0 


29 



Table 1. differential cross section F = £1X2 [ ^^Jidx^ — ^^dx^dx^ ^ ^ or a cnarm 
quark at large Ay > 4.0 and different ^/s. 



()2.24j) with pure LO DGLAP ones as a function of yfs at large rapidity gap 
Ay = 4.0. In the high energy region ^fs £ (1800GeV; 14TeV) LO DGLAP 
+ ln 2 x q va i give larger value of the cross section j^[c^(PP) — cr(PP)] than 
the pure LO DGLAP q va \ approach. This results simply from larger values 
of q va i distribution functions in LO DGLAP + ln 2 x approach at small x 
(x < 10- 2 ) in comparison to those obtained via LO DGLAP (or even NLO) 
evolution method. However double logarithmic terms ln 2 x play a significant 
role not only in the very small x region. Indeed, for large Ay one deals with 
not too small x in the integrals (|2.24|) or (|2.27|) because the lower limit for 
x is 

Xlow = ^Le^y (4.4) 
s 

In a case of a very large rapidity gap e.g. Ay = 8.0 at Tevatron energy 
y/s = 1800GeV and for = m c one should know parton distributions at 
x > 10 -2 instead of x > 3 • 10 -6 (Ay = 0). It could be very convenient 
because the very small x region is still experimentally unattainable. Large 
rapidity gap Ay implies however rapid decrease of suitable cross sections. 
So one should minimize s which is a suppression factor (s _1 ) keeping still 
large enough Ay ~ Ins. Our plots in Figs. 6-7 exhibit that for larger Ay 
the differential cross section [a(PP) — a(PP)] can change by almost 2 

(3) orders of magnitude from 81 nb (8 nb) at ^/s = 630GeV and Ay = 2.0 
to 1.2 nb (10 pb) at Ay = 7.0 and the same energy in a charm (bottom) 
production. For one value of Ay = 4.0 we can find j^[a(PP)—a(PP)] for a 
charm (bottom) case in the range from 22 nb (28 pb) at -y/i = lOOGeV to 2.8 
nb (340 pb) at y/s = 14TeV. Table I shows values of double differential cross 



15 



section x\X2[ — rfxirf-^, — -] for a charm quark at large Ay > 4.0 

and different yfs. These DGLAP + ln 2 x predictions which are equal from 
29 pb (Ay = 8.0, y/s = 14TeV) to 3.8 nb (Ay = 4.0, y/s = 39GeV) as not 
too small should probably be measurable. 



5. Summary. 

Difference of two cross sections a(PP) — a(PP)] for heavy quark pro- 
duction is in the Regge theory represented by pure mesonic reggeon ex- 
change process. In this way only valence quarks and antiquarks are taken 
into account in this hadron-hadron collisions. Distribution functions of va- 
lence quarks q va i(x, Q 2 ) which via factorization theorem enter to a suitable 
formula for quark-antiquark elastic scattering in cross sections can be ob- 
tained within PQCD approach. In our paper we have used a perturba- 
tive method which incorporates LO DGLAP evolution and resummation 
of double logarithmic terms ln 2 x as well. Using the dynamical GRV in- 
put parametrization at k 2 , = lGeV 2 we have found numerically in our LO 
DGLAP + ln 2 x approach q va i(x,Q 2 ) at higher scale Q 2 ~ 4m 2 and hence 
the cross section a(PP — > hh) — a(PP — > hh){s) and the differential cross 
section j^[<^(PP ~> hh) — a(PP — » hh)] for the charm and bottom produc- 
tion. Very interesting is an analysis where the rapidity gap Ay between the 
quark and the antiquark is large. The large rapidity gap (Ay > 2.0) causes 
that double logarithmic terms ln 2 x ~ Ay 2 become important. Processes 
with large rapidity gaps in high energy PP and PP collisions, observed 
at the Tevatron ^2] have been theoretically investigated e.g. in ^3], |14j . 
[T5| . [TT)] . The main reason that the large rapidity gap processes are so 
intensively studied is the opportunity to search for new particles e.g. Higgs 
bosons. Besides, the hard partonic processes with a rapidity gap between 
two jets in a final state enable examination of BFKL Pomeron behaviour 
|17j . |18j . Even though the pure mesonic exchange process gives a small 
contribution to cross sections in high energy hadron - hadron collisions, it 
could be interesting because it enable the PQCD test in double logarith- 
mic ln 2 x approximation. In our paper we have modified the mesonic Regge 
picture by PQCD LO DGLAP method combined with ln 2 x terms resum- 
mation. We have found that for the large rapidity gap Ay between a quark 
and an antiquark double logarithmic terms ln 2 x become significant. It must 
be however emphasized that in an unpolarized case it is very difficult to ob- 
serve a dominant role of the ln 2 x contribution to the flavour nonsinglet part 
of heavy quark production cross sections. There are two reasons of this fact. 
First: ln 2 x terms in q va l(x,Q 2 ), which by definition are significant in the 
small x region (x < 10 -2 ) are partially hidden behind leading Regge x a2 (°) 
behaviour of q va i- Second: the quark - antiquark scattering cross section 
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da/di is suppressed by a factor of s~ 2 , so one should minimize s with a 
reasonably large rapidity gap Ay = ln(s/Am 2 ). Larger Ay causes rapid 
decrease of the suitable cross sections, which however can be measured at 
Tevatron and LHC. 
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